Pole expansion of self-energy and interaction effect on topological insulators 
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We study eflFect of interactions on time-reversal-invariant topological insulators. Their topological 
indices are expressed by interacting Green's functions. Under the local self-energy approximation, 
we connect topological index and surface states of an interacting system to an auxiliary noninteract- 
mg system, whose Hamiltonian is related to the pole-expansions of the local self-energy. This finding 
greatly simplifies the calculation of interacting topological indices and gives an noninteracting pic- 
torial description of interaction driven topological phase transitions. Our results also bridge studies 
of the correlated topological insulating materials with the practical dynamical-mean-field-theory 
calculations. 
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Introduction-Recently, there has been a rising interest 
on a class of topological phase of matter: the topological 
insulators (TI), see [HE] for reviews. Questions like what 
is the interaction effect on the TI and how to characterize 
the topological phase beyond single particle basis pose in- 
triguing challenges. There are some previous studies on 
the interacting effect on TIs. Wang et al. |3] performed 
exact diagonalization and quantum Monte Carlo stud- 
ies of the interacting Haldane model. Hohenadler et al. 

Zheng et al. 0, Yu et al. |6i, Wu et al. i7j performed 
quantum Monte Carlo, variational cluster approximation 
and cluster dynamical-mean-field-theory (DMFT) stud- 
ies on the Kane-Mele-Hubbard (KMH) model. One in- 
teresting finding in these studies is that the topological 
phases could give their way to a gapped featureless state 
(spin liquid phase for KMH model) at small value of the 
spin orbital coupling. It is possible that the fluctuation 
effect driven by interaction breaks the topological phase 
significantly. This scenario goes beyond the mean-field 
treatments [3l |8]. Since in the latter case, long range 
orders (LRO) are usually needed to compete with the 
topological phases. 

However, one unsatisfactory feature in these studies is 
that there lacks a direct characterization of the topolog- 
ical order: the phase boundaries are mostly determined 
by indirect signatures such as the correlation functions 
for LRO or the closing and reopening features of excita- 
tion gaps, rather than a sharp change of one topological 
index. This leaves the mechanism and the nature of the 
correlation induced topological transitions unclear. 

In this Letter, we study the interaction effects on TI 
based on their single-particle Green's functions, we con- 
nect topological index and surface states of an interacting 
system to the indices of an auxiliary noninteracting sys- 
tem, whose Hamiltonian is related to the pole-expansions 
of the local self-energy. This finding greatly simplifies the 
calculation of interacting topological indices and gives an 
noninteracting pictorial description of interaction driven 



topological phase transitions. 

The system's topological index no is determined by 
the interacting Green's function through P| I10|: 

Gdf.G-'^Gd^G-^GdpG-'^GdaG-'^GdrG-^ (1) 

where G is Matsubara Green's function [22,, s^^p^T is 
the five-order anti-symmetric tensor, ^,v...t indices de- 
note the frequency- momentums {LL),kg;,ky,kz,kx). The 
same formula applies to the Chern number and Z2 in- 
dex. In latter case, fc^ is a pumping parameter extends 
three dimensional TI to four dimension[10 . Physically, 
Chern number describes quantum Hall conductance, Z2 
index describes magnetoelectric response of an insula- 
tor [T]. The insulator's Z2 index should be ether 
or 1 if there is time-reversal symmetry. In that case, 
only the parity of is meaningful. Formula Eqn|l] in- 
volves only the single particle Green's function and al- 
lows analysis of the effect of interaction, disorder and 
finite temperature on equal footing. In [11 , the au- 
thors point out the frequency-domain-winding-number 
(FDWN) could change the topological index without de- 
veloping long range orders. However, the paper made one 
crude assumption (besides the local self-energy assump- 
tions) that the self-energy is diagonal in the orbital basis 
and has no orbital dependence. In this paper, we relax 
this assumption by adopting the pole-expansion (PE) for 
general local self-energies. We link the topological in- 
dices and surface states of an interacting system to an 
noninteracting system thereof. 

Pole expansion of local self-energy- Considering a TI 
with noninteracting Hamiltonian H^, who is a matrix of 
the size m. Assume that the local interaction generates 
a self-energy without momentum dependence except in 
the Hartree-Fock part, thus, the self-energy can be writ- 
ten into the PE form as S(w, k) = Sk + V^{iuj - P)-^V. 
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This type of self-energy can be obtained by the DMFT. 
which is believed to be reliable for 3D correlated sys- 
tems [12] ■ Sk is a TO X m matrix denotes the asymptotic 
(Hartree-Fock) contributions to the self-energy. P and V 
are frequency-independent matrices of the size Np x Np 
and Np x to respectively. P and V determine the po- 
sition and weight of poles of the self-energy. Physically, 



PE corresponds to represent the interaction induced local 
self-energy to coupling with noninteracting baths. One 
could introduce a noninteracting pseudo-Hamiltonian of 

the form ^ Hu = (^^^ ^ ^ ^ ^p^, which has the 

matrix size N = Np + to. The corresponding Green's 
function reads 



G = {iuj- Hk)-^ = 



-V iuj-P 



[iiu + /I - i/k - Sk - Vfiiiu - P)-^V]-^ 



Clearly, G is related to the interacting Green's function G = [iuj+fi — H]f^ — J^{oj)] ^ through a projection G = V^GV. 

/I ., 
1 ., 



The projection operator V = 
the integrand of Eqn{T]is: 



is a X TO matrix projects to the upper-left elements of G. Notice that 



V 



TrGd^G- ^Gdk^G-^Gdk^G-^Gdk^G-^ Gdk, ^ 
BG 

= -TY—dk^G~^Gdk^G-^Gdkfi-^Gdk^G-^ 
BG 

= -TiV'^—VV^dk^G-^VV^GVV^dk^G-^VV'^GVV'^BkG-^VV'^GVV^dk^G-^V 

OUJ 

dC 

= -Tr—VV''Bk^G-^VV^GVV^Bk,G-^VV^GVV^Bk^G-'^VV^GVV^dk^G-'^VV^ 

BU! " 

BG 

= -TT—Bk^G-'Gdk^G-^GBk^G~^Gdk,G-^ 

= TvGB^G-^Gdk^G'^GBk^G-^GBk^G'^GBk.G-' (2) 



where we have used the identities dk^G ^ = — 9fei(i?k + 
Sk) = V^dk,G-^V a.iidVV^ dk^G-^Vr^ = Bk^G-^- Base 
on Eqnj2] one can proof that the topological (Chern or 
Z2) indices calculated from G and G are identical, i.e., 

no = (3) 

Further more, the surface and bulk Green's functions 
are related through the projection G^™^ = Q^GQ. Q 
consists of a Fourier transformation of (let's say) momen- 
tum kz to real space, then project to the surface norms 
to 2;-axis. Since the two projection operators V and Q 
act in different spaces, one has ['P,Q] =0. Thus, 

g,surf ^ Qt'ptG-pQ = r^Q^GQP = ■P^G'^^T' (4) 

which means that the surface Green's function of an 
interacting system {G^^"^) is related to an auxiliary 
noninteracting system's surface Green's function (G^"''^) 
through the same projection V. 



Eqnj3]and|4]are the central results of the paper. They 
link the topological index and surface states of an in- 
teracting model to a noninteracting model. The lat- 
ter is purely determined by the noninteracting pseudo- 
Hamiltonian i^k- For a given local self-energy, one could 
construct its corresponding pseudo-Hamiltonian and cal- 
culate its topological index fellowing conventional pro- 
cedures for noninteracting systems [T4l - fT6| . This fact 
greatly reduces the computational complexity of calcu- 
lating topological indices of interacting systems. More- 
over, the relationship Eqnj3]and Eqn|4] provide a way of 
understanding interaction induced topological transition 
by its noninteracting counterpart. 

Applications- We demonstrate how the general rules 
Eqn|3]and Eqn|4] applies to an concrete model for inter- 
acting Z2 insulator. We choose noninteracting Hamilto- 
nian to be |23| . 

ffk = sinfc^r^ + sin kyV^ + sin/c^r^ + 7W(k)r^ (5) 

where A^(k) = to — 3 -f- {cosk^ + cos ky + coskz). The 
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FIG. 1: (Color online) Z2 index (red line) and gap (blue line) 
of noninteracting model (Eqn[5| as a function of m. The phase 
boundaries between different Z2 indices are in accordance with 
the gap closing points. 

Hamiltonian describes a 3D (strong) Z2 TI for < m < 2 
and 4 < TO < 6. Inside the TI phase the gap of the system 
A scales linearly with m, see Figjl] We set to = 1 in the 
following studies. Notice that the eigenvalues of iJk are 

iii^k, with i?k = ^sin^ + sin^ ky + sin^ + M'^(k). 
It also determines the band gap A = 1 and band top 
r = 5. 

In the following we consider the effect of local interac- 
tion on the TI phase. Rather than trying to solve the in- 
teracting model explicitly (leave to future studies), we set 
up a few ansatz of the self-energies (see Tablejl| and cal- 
culate its interacting Z2 index using the following meth- 
ods: (1) We get the interacting Z2 number based on its 
FDWN [U]; (2) With the PE coefficients of the given self- 
energy (see Tabl^ we construct the pseudo-Hamiltonian 
Hk and calculate its Z2 index fellowing noninteracting 
procedures jTHj. We also study the gap closing condition 
and surface states of the interacting system through the 
noninteracting reference system H^. 

We start with the case (a): T,(uj, k) = ji^zrp- We draw 

the path defined by cj i— > G^^lp{uj) = iuj — j^irp on the 
complex plane and find that the topological transition 
point corresponds to where the path touch the lower edge 
of the eigenenergy continuum, i.e. the band gap A. This 
defines the phase boundary between the TI and metal 
phase ^ = ±A, see Figj2](a). We then calculate the Z2 
index of the noninteracting pseudo-Hamiltonian iJk with 
its single-particle engenstates [23] • The resulting phase 
diagram is shown in Figj2jb). Phase boundaries are in 
accordance with the prediction based on FDWN. 

The surface Green's function could be directly calcu- 
lated for the interacting system with given self-energy 
or from the surface Green's function associated to the 
pseudo-Hamiltonian H-^, (Eqn|4]). We have numerically 
verified that the two approaches give identical results. 
The resulting surface spectral functions are shown in 
Figj3] There are gapless surface states crossing the Fermi 
surface when the system is in the interacting topological 
phase. Fig |3] upper panel. 

We then move to the case (b): Sfa;, k) = . ^, „ + . ^ „ 
|25j . It has been shown that the interacting system is 
gapful and its Z2 index equals to for P = and 1 




FIG. 2: (Color online) (Left panel) Paths of the mapping 
^ ^ G^r^p ^ ii^ - with y = 1, P = 2 (blue line) 

and P = 0.15 (red line). Gray line indicates the large circle 
we appended to close the contour. Crosses on real axis indi- 
cate the eigenenergy continuum of the noninteracting Hamil- 
tonian, whose lower (upper) bound equals to A (F). Topolog- 
ical transition between TI and metal occurs when |P| A — , 
where the blue line touches the lower bound. The gap reopens 
for |P|r < , results into a normal insulating phase. (Right 
panel) Phase diagram determined by calculating the Z2 index 
of the pseudo-Hamiltonian (see the PE coefficients in Table|l| . 
Shaded and white region have Z2 = 1 and respectively, the 
yellow region denotes the metallic state. The phase bound- 
aries |P|A — and |P|r — are in accordance with the 
FDWN analysis. 

otherwise [H]. Diagonalizing i/k shows that the nonin- 
teracting system is gapful for P ^ 0, with Z2 index equal 
to 1, which validates the aforementioned result. The gap 
of the pseudo-Hamiltonian closes at P = 0. At first look, 
this seems to be inconsistent with the Eqnj3] because the 
ground state of is gapless while the interacting sys- 
tem is gapful. This paradox is solved by noticing that for 
P = the two poles of case (b) merge into one but we are 
still constructing the pseudo-Hamiltonian with Np = 8. 
When treat this particular point with Np = 4, case (b) 
reduced to case (a) with P = 0, which is topological triv- 
ial and consistent with the FDWN prediction. Lesson 
from this example is that when constructing the pseudo- 
Hamiltonian, degenerate poles should be treated care- 
fully. 

Discussions- Direct dealing with the frequency- 
momentum integration in Eqn[T]is the most general way 
of calculating interacting topological index. It does not 
make any priori assumptions on the form of the self- 
energies. However, as been pointed out in |11| . the distri- 
bution of the integrand in Eqn^is very inhomogeneous. 
Numerical integration for inversion symmetric Z2 TIs is 
even harder, since all of the topological information con- 
centrate around those time-reversal invariant momenta 
|17| . Integrate with Monte Carlo technique or adap- 
tive meshes may overcome the difficulties to some de- 
gree. On the other hand, the concept of FDWN greatly 
simplify the calculation and gives an appealing picture 
for the interaction effect on TIs. However, the useful- 
ness of FDWN is only limited to the case where the self- 
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TABLE I: Self-energies and their PE coefBcients. Interacting Z2 indices calculated by two methods (see text) are in accordance 
with each other. 



Cases S(a; k) Np Ek P Z2 = 1 Z2 = 

~{^) 4 vh, PI4, |P|A > |P|r < y2 

(b) T^ + i^ 8 {V,V)(g)U diag{P,-P)(E)U P/0 P = 




FIG. 3: Surface spectral functions for self-energy case (a) with 
V = 1. (upper panel) P = 2 is in TI phase, (lower panel) 
P — 0.15 is normal insulator. The chemical potential lies at 
zero. 



For two dimensional systems, ground state wave func- 
tion (GSWF) under twisted boundaries could be used to 
extract topological index[TB]. However, GSWF is very 
hard to get for higher dimensional interacting fermionic 
models. In fact, comparing to Green's functions, GSWF 
contains redundant information for topological proper- 
ties. Furthermore, there are large number of analytical 
or numerical approaches of getting interacting Green's 
function for correlated fermionic models. 

Base on these analysis, we anticipate the local self- 
energy approximation and PE technique would of great 
use in subsequent studies of 3D interacting TIs. 

Conclusion- We show that the local self-energy con- 
tains the topological signature of interacting TIs. With 
the help of pole-expansion of the local self-energy, we 
find that the topological index of an interacting topolog- 
ical (Chern or Z2) insulator is identical to the one for an 
noninteracting model. This findings greatly simplify the 
calculation of the topological index for interacting TIs 
and give an appealing noninteracting picture for interac- 
tion induced topological phase transitions. Techniques 
report here could have greatly potential use in direct 
determination of the topological indices in the realistic 
LDA+DMFT studies of the 3D correlated TI materials. 

Note- See Ref.[T^ for a recent application of pole- 
expansion method on Bernevig-Hughes-Zhang model. 
There is an independent application of FDWN and pole- 
expansion idea to topological quantum phase transitions 
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energy is diagonal in the orbital basis (besides the local 
assumptions). The PE approach lies in the middle of the 
two aforementioned limiting cases. It relaxes the orbital- 
independent assumption and is more general for higher 
dimensional TIs. PE is merely the discrete approxima- 
tion for the self-energy, which provides a systematic way 
to approach the true self-energy by making the number 
of poles denser and denser. 

Under the local self-energy approximation, the PE ap- 
proach could also be used to calculate the Chern and 
Z2 indices for two dimensional systems. However, the 
main advantage of PE approach lies in high dimensions. 
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